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We report on a search for pair production of first-generation scalar leptoquarks (LQ) in pp colli- 
sions at -^=1.96 TeV using an integrated luminosity of 252 pb _1 collected at the Fermilab Tevatron 
collider by the D0 detector. We observe no evidence for LQ production in the topologies arising 
from LQLQ — » eqeq and LQLQ — > eqvq, and derive 95% C.L. lower limits on the LQ mass as a 
function of f3, where (3 is the branching fraction for LQ — > eq. The limits are 241 and 218 GeV/c 2 
for (3=1 and 0.5, respectively. These results are combined with those obtained by D0 at v / s=1.8 
TeV, which increases these LQ mass limits to 256 and 234 GeV/c 2 . 

PACS numbers: 14.80.-j, 13.85.Rm 



Several extensions of the standard model (SM) include 
leptoquarks (LQ) which carry color, fractional electric 
charge, and both lepton (I) and quark (q) quantum num- 
bers and would decay into a lepton and a quark . The 



HI and ZEUS experiments at the e^p collider HERA at 
DESY published Q lower limits on the mass of a first- 
generation LQ that depend on the unknown leptoquark- 
l-q Yukawa coupling A. At the CERN LEP collider, pair 
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production of leptoquarks could occur in e+e~ collisions 
via a virtual 7 or Z boson in the s-channel. At the Fer- 
milab Tevatron collider, leptoquarks would be pair pro- 
duced dominantly through qq annihilation (for Mlq > 
100 GeV/c 2 ) and gluon fusion. Such pair production 
mechanisms are independent of the coupling A. Experi- 
ments at the LEP collider and at the Fermilab Teva- 
tron collider [3 00 set lower limits on the masses of 
leptoquarks. In this Letter, we present a search for first- 
generation scalar leptoquark pairs produced in pp colli- 
sions at -^5=1. 96 TeV for two cases: when both lepto- 
quarks decay to an electron and a quark with a branch- 
ing fraction (Br) (3 2 , where (3 is the leptoquark branching 
fraction into an electron and a quark, and when one of the 
leptoquarks decays to an electron and a quark and the 
other to a neutrino and a quark with Br = 2/3(1-/3). The 
final states consist of two electrons and two jets (eejj) 
or of an electron, two jets, and missing transverse energy 
corresponding to the neutrino which escapes detection 
(evjj)- 

The D0 detector y| comprises three main elements. 
A magnetic central-tracking system, which consists of a 
silicon microstrip tracker and a central fiber tracker, is 
located within a 2 T superconducting solenoidal magnet. 
Three liquid-argon/uranium calorimeters, a central sec- 
tion (CC) covering pseudorapidities r\ 8] with \rj\ up to 
« 1 and two end calorimeters (EC) extending coverage to 
\rj\ rj 4 jjj, are housed in separate cryostats. Scintillators 
between the CC and EC cryostats provide sampling of 
developing showers for 1.1 < \rj\ < 1.4. A muon system 
is located outside the calorimeters. 

The data used in this analysis were collected from April 
2002 to March 2004. The integrated luminosity for this 
data sample is 252 ± 16 pb _1 . Events were required to 
pass at least one of a set of electron triggers based on 
the requirement of one electromagnetic trigger tower to 
be above threshold and on shower shape conditions. The 
efficiencies of the trigger combinations used in the eejj 
and evjj analyses have been measured using data. They 
are ~ 100% for two electrons of transverse energy (E^ M ) 
above 25 GeV, and for one electron above 40 GeV. The 
small loss of events due to the trigger inefficiencies for 
E® M below 40 GeV is taken into account using proper 
weighting for Monte Carlo (MC) events. 

Electrons are reconstructed as calorimeter electromag- 
netic clusters which match a track in the central-tracking 
system. Electromagnetic (EM) clusters are identified 
by the characteristics of their energy deposition in the 
calorimeter. Cuts are applied on the fraction of the en- 
ergy in the electromagnetic calorimeter and the isolation 
of the cluster in the calorimeter. EM clusters are marked 
as tight when they satisfy a shower shape condition and 
loose otherwise. Jets are reconstructed using the iter- 
ative, midpoint cone algorithm jicij | with a cone size of 
0.5. The energy measurement of the jets has been cali- 
brated as a function of the jet transverse energy and 77 by 



balancing energy in photon plus jet events. The missing 
transverse energy (fir) is calculated as the vector sum of 
the transverse energies in the calorimeter cells, removing 
contributions from detector noise. 

For both channels, the background arising from mul- 
tijet events is determined from a sample of data events 
(QCD sample) that satisfy the main cuts used in the 
analysis except that each EM cluster is loose instead of 
tight. A QCD normalization factor is extracted for this 
sample in a part of the phase space where the LQ con- 
tribution is expected to be negligible. The QCD sample 
normalized by this factor is used to derive the multijct 
contribution in the relevant part of the phase space. To 
evaluate the Z boson/Drell-Yan (Z/DY) and the W bo- 
son background contributions, samples of MC events gen- 
erated with ALPGEN [ll| or PYTHIA were used. Sam- 
ples of pythia tt events (m t = 175 GeV/c 2 ) were used 
to calculate the top quark background. LQLQ — ► eejj 
and LQLQ — ► evjj MC samples were generated using 
PYTHIA for LQ masses from 120 to 280 GeV/c 2 in steps 
of 20 GeV/c 2 . All MC events were processed using a 
full simulation of the detector based on GEANT and 
the complete event reconstruction. The efficiencies of the 
various cuts, measured using the data, were taken into 
account using proper weightings of the MC events. 

The eejj analysis requires two tight EM clusters with 
Ef^ > 25 GeV and at least two jets with E T > 20 GeV 
within \rj\ < 2.4. At least one of the EM clusters should 
spatially match an isolated track and at least one should 
be in the CC fiducial region. The major SM background 
sources that mimic the eejj decay of a LQ pair are multi- 
jet events (where two of the jets are misidentified as EM 
objects), Z/DY production, and top quark pair produc- 
tion. To suppress background from Z boson production, 
events with a di-electron mass (M2EM.) compatible with 
the Z boson mass (80 GeV/c 2 < M 2EM < 102 GeV/c 2 ) 
are rejected. Finally St > 450 GeV is also required, 
where St is the scalar sum of the transverse energies of 
the two electrons and the two leading jets. In Fig.Qt, the 
St distributions for data and background after applying 
the Z boson mass cut are shown. This choice of the cut- 
off has been optimized using MC signal and background 
events to get the best expected mass limit. The total 
efficiencies for a LQ signal are summarized in Table [I] 
The multijet background is estimated from two samples 
of events with two EM clusters M > 15 GeV which 
have at least one matched track and no reconstructed 
jets. Both EM clusters are tight in one sample and loose 
in the other. The QCD normalization factor is deter- 
mined by the normalization of the M2EM distributions of 
the two samples below 75 GeV/c 2 . The Z/DY and top 
quark contributions are normalized to the integrated lu- 
minosity. Table ^ lists the number of events in the data 
and the number of expected events from SM background 
sources. 

Systematic uncertainties on the background are deter- 
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TABLE I: Efficiencies after ali cuts and 95% C.L. upper limits 
on production cross section x branching fraction Br, as a 
function of A/lq, for the two channels. 



M LQ (GeV/c 2 ) 


eejj 




evjj 






e(%) axBr(pb) 


e(%) axBr(pb) 


120 


2.2±0.5 


0.950 


4.6±0.5 


0.34 


140 


4.5±0.9 


0.444 


7.9±0.8 


0.20 


160 


8.9±1.7 


0.223 


11.7±1.1 


0.14 


180 


12.6±2.4 


0.156 


15.5±1.5 


0.10 


200 


18.5±3.0 


0.102 


17.8±1.7 


0.088 


220 


24.6±3.5 


0.075 


18.9±1.8 


0.083 


240 


30.3±3.9 


0.060 


20.9±1.9 


0.075 


260 


34.0±4.0 


0.053 


21.9±2.1 


0.071 


280 


36.0±4.0 


0.050 


22.7±2.1 


0.069 
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FIG. 1: The St distributions for the eejj events (a) and 
evjj events (b) from data (triangles) compared to the SM 
background (solid histograms). The dot-dashed histograms 
are the expected distributions for a 240 GeV/c 2 LQ signal (a) 
and for a 200 GeV/c 2 LQ signal (b). 



mined to be 15% from the QCD normalization factor and 
6% from the efficiencies of the identification of electrons 
and jets (particle- ID). An uncertainty (26%) from the 
jet energy scale is determined by varying the correction 
factor on the calorimeter response to jets by one stan- 
dard deviation. A systematic uncertainty on the Z /DY 
background (20%) is calculated by taking into account 
the differences between the two Z /DY MC samples. On 
the signal, the particle-ID and the limited statistics of 
the MC sample correspond to systematic uncertainties 
of 6% and 1.2% respectively. Comparing acceptances for 
the signal samples generated with pythia using different 
parametrizations of parton distribution functions (PDFs) 
leads to an uncertainty of 5%. The uncertainty due to 
the jet energy scale is dependent on the LQ mass (7.3% 
for a LQ mass of 240 GeV/c 2 ). The total uncertainty 
on the efficiency is (17-9)% in the mass range 180-280 
GeV/c 2 . 

The data are consistent with the expected SM back- 
ground and no evidence for leptoquark production is ob- 
served in the eejj channel. Thus we can set an upper 
limit at the 95% C.L. on the LQ pair production cross 
section using a Bayesian approach |b|. The limits are 
tabulated in Table |U and shown in Fig. [2^ as a function 
of LQ mass. To compare our experimental results with 
theory, we use the next-to-leading order (NLO) cross sec- 



TABLE II: Number of events in data compared with back- 
ground expectation at different stages of the eejj analysis. 



eejj 



Z boson veto S T >450 GeV 



Data 
Total background 
Z/DY + jets 
_ Multijet 
tt production 



467 
406 ± 100 
342 ± 99 

59 ±16 
4.7 ±0.4 



95 
92±17 
41 ± 11 
47 ± 13 
3.8 ±0.3 



1 

0.54±0.11 
0.22±0.07 
0.27±0.08 
0.05±0.01 




180 200 220 240 260 280 

Scalar Leptoquark Mass (GeV/c 2 ) 




120 140 160 180 200 220 240 260 280 

Scalar Leptoquark Mass (GeV/c 2 ) 



FIG. 2: The 95% C.L. limit on the experimental cross section 
times branching fraction as a function of LQ mass (circles) 
for the eejj (a) and evjj (b) channels. The NLO theoreti- 
cal cross sections [l5| are plotted for different values of the 
renormalization scale factor: Mlq (full line), Mlq/2 (dot- 
ted curve) and 2Mlq (dashed curve) taking into account the 
PDF uncertainties. A mass limit of 241 GeV/c 2 (a) and of 208 
GeV/c 2 (b) for first-generation scalar leptoquarks is obtained 
for P=l and /3=0.5, respectively. 



tion for scalar leptoquark pair production from Ref . [15| , 
with the CTEQ6 PDF [ljj. The theoretical uncertainties 
correspond to the variation from Mlq/2 to 2Mlq of the 
renormalization scale /i used in the calculation and to the 
errors on the PDFs. To set a limit on the LQ mass we 
compare our experimental limit to the theoretical cross 
section for /i = 2Mlq, which is conservative as it corre- 
sponds to the lower value of the theoretical cross section. 
The value of the theoretical cross section would increase 
by ~ 7% if the PDF errors were neglected. A lower limit 
on the leptoquark mass of 241 GeV/c 2 is obtained for 
(3=1. 

The evjj analysis requires exactly one tight EM clus- 
ter (£f M > 35 GeV) in the CC fiducial region which 
matches an isolated track spatially and kinematically. At 
least two jets with Et > 25 GeV within I77I < 2.4 and 
|?t > 30 GeV are required. The main SM background 
sources which would mimic the evjj decay of a LQ pair 
are events with multijet production (where a jet is re- 
constructed as an electron and the Ifir comes from jet 
mismeasurements) , W + 2 jets events, and top quark 
pair production. A veto on muons with px > 10 GeV/c 
is applied to reduce the di-lepton background from tt 
decays. To suppress background from W boson produc- 
tion, events with a transverse mass of the electron and 
the missing energy M T V < 130 GeV/c 2 are rejected. Fi- 
nally St > 330 GeV is required, where here St is the 
sum of the transverse energies of the electron, the two 
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TABLE III: Number of events in data compared with back- 
ground expectation at different stages of the evjj analysis. 
The values of the cuts are in GeV or in GeV/c 2 . 





$t > 30 


Mr > 130 


St 


> 


330 


Data 


900 


14 




1 




Total background 


902 ±211 


13.9 ±4.4 


3.6 


± 


1.2 


W + jets 


811 ±211 


10.0 ±4.4 


2.2 


± 


1.2 


Multijet 


76 ±7 


2.3±0.5 


0.72 


± 


0.28 


tt production 


14.7±2.9 


1.6 ±0.37 


0.70 


± 


0.17 



jets, and the ftp. The distribution of the variable St for 
the data and the total background is shown in Fig. 2J) 
after applying the M T V cut. The choice of the cutoff has 
been optimized as above. The total efficiency of these 
cuts for a LQ signal is given in Table |U To determine the 
multijet background we use a data sample that passed all 
the preceeding cuts but with a loose EM cluster match- 
ing spatially a track. The QCD normalization factor is 
determined using the ratio of the number of events with 
fx < 10 GeV in this and in the search samples. The W 
boson background is normalized to the data at transverse 
mass 60 GeV/c 2 < Mf < 100 GeV/c 2 . The top quark 
background is normalized to the integrated luminosity 
using the NNLO theoretical cross section. The number 
of events which survive the cuts and the number of pre- 
dicted background events are summarized in Table InTl 

Systematic uncertainties associated with the QCD nor- 
malization factor (9%) and W boson normalization fac- 
tor (5.7%) are determined by the limited statistics of the 
samples and the choice of kinematical domain over which 
the normalization is done. The jet energy scale uncer- 
tainty introduces uncertainties equal to 25% for W boson 
production and 8.5% for the top-quark-pair production. 
For the W boson background an uncertainty equal to 33% 
is associated to the shape of the $t distribution. A 25% 
error has been included as systematic uncertainty on the 
top quark cross section. Finally, there is an uncertainty 
of 3.8% on the particle-ID acceptance. Three systematic 
uncertainties are determined on the signal acceptance: 
3.8% comes from the uncertainty on the particle-ID, 5% 
is due to the jet energy scale uncertainty, and 5.4% cor- 
responds to the acceptance variations for different PDF 
paramcterizations. 

As no excess of data over background is found in the 
evjj channel, an upper limit on the production cross sec- 
tion for a first-generation scalar leptoquark is derived and 
shown in Fig. and in Table [i] A comparison of these 
limits to theoretical calculations of the cross section , 
performed as described above, gives a lower limit on 
the first-generation scalar LQ mass of 208 GeV/c 2 for 
(3 = 0.5. 

Combination of the limits obtained in the searches in 
the eejj and evjj channels is done using a Bayesian like- 
lihood technique |l7j , with correlated uncertainties taken 



TABLE IV: 95% C.L. lower limits on the first-generation 
scalar leptoquark mass (in GeV/c 2 ), as a function of {3. The 
mass limits from D0 {eejj, evjj and vvjj combined) |4| and 
CDF {eejj) || at Run I (~ 120 pb _1 ) are also given, as well 
as the limits obtained by combining the D0 Run I and Run 
II results. 

P 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1. 

eejj 158 180 203 220 232 241 

evjj 169 193 203 207 208 207 203 193 169 

Comb. Run II 169 193 204 212 218 223 228 232 237 241 
D0 Run I 110 204 225 

D0 Runs I & II 183 206 218 227 234 239 244 248 252 256 
CDF Run I 213 




140 160 180 200 220 240 260 



Scalar Leptoquark Mass (GeV/c ) 

FIG. 3: Excluded regions (shaded area) at the 95% C.L. in the 
f3 versus LQ mass plane for the production of first-generation 
scalar leptoquarks. 

into account. The limits on the cross sections obtained 
at the 95% C.L. for the combination of the two channels 
and different values of (3 are compared with the NLO LQ 
pair production cross section [lFj and lower mass limits 
are derived and given, as a function of /3, in Table ITVl and 
shown in Fig. |3| In Table |IV| are also shown the Run I 
mass limits based on an integrated luminosity ~ 120 pb -1 
obtained by D0 Q, using the three channels eejj, evjj 
and vvjj, and CDF [j| {eejj channel). This analysis 
sets a 95% C.L. limit on the first-generation leptoquark 
mass of M L q > 218 GeV/c 2 for £=0.5, and M L q > 241 
GeV/c 2 for /3=1. The D0 Run II and Run I results are 
combined, using the same method, and the results are 
shown in Table |IV] and in Fig. The 95% C.L. limits 
on the first-generation leptoquark mass are Mlq > 234 
GeV/c 2 for /3=0.5, and Af LQ > 256 GeV/c 2 for /3=1. 
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